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Abstract 

Using our recently developed Chebyshev expansion technique for finite-temperature dynamical correlation functions we 
numerically study the AC conductivity aiuj) of the Anderson model on large cubic clusters of up to 100^ sites. Extending 
previous results we focus on the role of the boundary conditions and check the consistency of the DC limit, a; ^ 0, by 
comparing with direct conductance calculations based on a Greens function approach in a Landauer Biittiker type setup. 
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The numerical calculation of linear response functions is 
one of the standard tasks in condensed matter theory and 
many other areas of physics. In practice, however, the di- 
mension A'^ of the corresponding eigenvalue problem usually 
becomes enormously large. A complete diagonalisation and 
a naive evaluation of linear response functions is prohibitive 
is such situations, since the required time would scale at 
least as A^'^. In a recent work [1] one of us proposed an ad- 
vanced Chebyshev expansion method for the calculation of 
dynamical correlation functions at finite temperature, which 
is linear in the system size N . To yield data for all tem- 
peratures and chemical potentials within a single simulation 
run marks one of its additional advantages. As a particu- 
larly interesting application we considered the optical (AC) 
conductivity (j{lo) of non-interacting electrons in a random 
potential, which can be described by the well-known An- 
derson model [2] 
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Here c[ (c^) denote fermion creation (annihilation) operators 
and £ [~W/2, W/2] is the uniformly distributed local 
potential. 
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In the present work we extend our recent studies and 
focus on the sensitivity of a{uj) to changes in the boundary 
conditions. In addition, we check the consistency of the 
w — > limit by comparing extrapolations of the AC data 
and direct calculations of the DC conductance combined 
with finite size scaling. 

Within linear response theory the real part of the optical 
conductivity is given by 
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where |n) and |m) denote eigenstates of the Hamiltonian 
with energies E„ and Em, i^nm = En — Em, f{E) = 
[exp{/3{E — ft)) + is the Fermi function, and Jx — 
—if^^iiclci+x ^ '^i+x'^i) '^he a::-component of the current 
operator. Even at zero temperature Eq. (2) involves a sum- 
mation over matrix elements between all one-particle eigen- 
states of H, which for reasonably large systems can hardly 
be calculated explicitly. Consequently, until now, the num- 
ber of numerical attempts to this problem is very limited [3]. 
Rewriting (j{lu) in terms of Fermi integrals over a matrix 
element density j{x, y), which itself is independent of tem- 
perature and chemical potential, our new approach [1] is 
based on the observation that the latter can be efficiently 
calculated with a two-dimensional generalisation of Cheby- 
shev expansion and the Kernel Polynomial Method [4]. This 
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Fig. 1 . Optical conductivity of the 3D Anderson model at T — and ^ — 
(band centre) for increasing disorder W. Bold lines mark W/t — 16, 
which approximately conesponds to the critical disorder. Most curves 
coirespond to — 50"^ sites, expansion order M — 1024 and averaging 
over 240 samples; sets that extend to lower uj are based on A'^ — 100^, 
M = 2048 and 400 samples. 




Fig. 2. Finite size scaling of the DC conductivity calculated directly [5] 
within the Landauer Blittiker type geometry shown in the inset. 

leads to an 0{N) algorithm and allows the study of systems 
with as many as 100'^ sites. 

The good quality of the obtained conductivity data allows 
for tests of the various analytical predictions for the finite 
frequency behaviour. In particular, for d > 2 dimensions 
we expect Aa = o-{uj) — (t(0) ~ ti;('*~2)/2 metallic 
phase, cr(to) ~ u^''-^^)/'^ at the metal-insulator transition, 
and cr(ci;) '-^ o;'^ far inside the insulating phase [6]. Figures 1 
and 3 illustrate that for low frequency cr(cj) is indeed well 
described by a power law cr{0) + C to" . As expected, for 
samples with open boundary conditions (OBC) the extrap- 
olated DC contribution cr(0) is always zero, whereas for 
periodic (PBC) and anti-periodic (ABC) boundary condi- 
tions we observe a continuous metal-insulator transition at 
W/t « 16. As a further check we also calculated cr(0) di- 
rectly using a numerical Greens function approach within a 
Landauer Biittiker type geometry [5] and the finite size scal- 
ing ansatz cr(0) = Le'^'^ /[a + L) (cf. Fig. 2). Interest- 
ingly, within both approaches the resulting critical exponent 
for (t(0) seems to be of the order of 2, larger than the ex- 
pected V « 1.57, a puzzle which certainly requires further 
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Fig. 3. Exponent a and DC conductivity ct(0) obtained from fits of the 
low-frequency conductivity to <t{uj) — + Ccj" (vertical bars in 

Fig. 1 mark the underlying frequency range). En"or bars are estimated by 
slightly vaiying G [-0.051V, 0.05W]. 

examination. Focussing on the exponent a, our PBC and 
ABC data appears to be consistent with the analytic predic- 
tions, i.e., a 0.5 for weak disorder, a ~ 1/3 at the criti- 
cal disorder, and eventually a ^ 2 for very strong disorder. 
Note, however, that for OBC a is a simple, monotonously 
increasing function of W . The above results show, that the 
new technique is a reliable numerical tool for large scale 
calculations of dynamical correlations functions. 

We acknowledge financial support by the Australian Re- 
search Council, calculations were performed at APAC Can- 
berra, ac3 Sydney and NIC Jiilich. 



References 

[1] A. WeiBe, Eur Phys. J. B 40 (2004) 125. 
[2] P. W. Anderson, Phys. Rev. 109 (1958) 1492. 

[3] R. C. Albers, J. E. Gubernatis, Phys. Rev. B 17 (1978) 4487; A. Singh, 
W. L. McMillan, J. Phys. C 18 (1985) 2097; M. Hwang, A. Gonis, 

A. J. Freeman, Phys. Rev. B 35 (1987) 8974; P Lambrianides, H. B. 
Shore, Phys. Rev. B 50 (1994) 7268; H. Shima, T. Nakayama, Phys. 
Rev. B 60 (1999) 14066. 

[4] R. N. Silver, H. Roder, Int. J. Mod. Phys. C 5 (1994) 935. 

[5] J. A. Verges, Comp. Phys. Comm. 118 (1999) 71. 

[6] N. F Mott, Adv Phys. 16 (1967) 49; F J. Wegner, Z. Phys. B 25 
(1976) 327; R. Oppermann, F Wegner, Z. Phys. B 34 (1979) 327; 

B. Shapiro, E. Abrahams, Phys. Rev. B 24 (1981) 4889; B. Shapiro, 
Phys. Rev. B 25 (1982) 4266. 



2 



